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PLANET FORMATION IN BINARY SYSTEMS: A SEPARATION-DEPENDENT MECHANISM? 
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ABSTRACT 

In this article, I examine several observational trends regarding protoplanetary disks, debris disks 
and exoplanets in binary systems in an attempt to constrain the physical mechanisms of planet 
formation in such a context. Binaries wider than about 100 AU are indistinguishable from single stars 
in all aspects. Binaries in the 5-100 AU range, on the other hand, are associated with shorter-lived 
but (at least in some cases) equally massive disks. Furthermore, they form planetesimals and mature 
planetary systems at a similar rate as wider binaries and single stars, albeit with the peculiarity 
that they predominantly produce high-mass planets. I posit that the location of a stellar companion 
influences the relative importance of the core accretion and disk fragmentation planet formation 
processes, with the latter mechanism being predominant in binaries tighter than 100 AU. 
Subject headings: binaries: general — planetary systems — planetary systems: formation — planetary 
systems: protoplanetary disks 



1. INTRODUCTION 

The exponentially growing number of known extraso- 
lar planets now enables statistical analyses to probe their 
formation mechanism. Two theoretical frameworks have 
been proposed to account for the formation of gas gi- 
ant planets: the slow and gradual core accretion model 
(Lissauer & Stevenson 2007), and the fast and abrupt 
disk fragmentation model (Durisen et al. 2007) . The de- 
bate regarding their relative importance is still ongoing. 
Both mechanisms may contribute to planet formation, 
depending on the initial conditions in any given pro- 
toplanetary disk (Boley 2009, and references therein). 
By and large, our understanding of the planet formation 
process is focused on the case of a single star+disk sys- 
tem. Yet, roughly half of all solar-type field stars, and 
an even higher proportion of pre-main sequence (PMS) 
stars, possess a stellar companion (e.g., Duquennoy & 
Mayor 1991; Mathieu 1994; Duchene et al. 2007). Since 
the disk and multiplicity phenomena are associated with 
similar ranges of distances from the central star, the dy- 
namical influence of a companion on a disk may be dra- 
matic. Theory and observations agree that stellar com- 
panions can open large gaps in disks, or truncate them 
to much smaller radii than they would otherwise have 
(e.g., Artymowicz & Lubow 1994; Ireland & Kraus 2008). 
The consequences for planet formation are still uncertain, 
however. 

Observations of protoplanetary disks among PMS stars 
have revealed that tight binaries generally show substan- 
tially reduced (sub) millimeter thermal emission (Beck- 
with et al. 1990; Jensen et al. 1996) as well as a much 
rarer presence of small dust grains in regions a few AU 
from either component (Chez et al. 1993; Cieza et al. 
2009). Both trends can be qualitatively accounted for 
by companion-induced disk truncation, which can simul- 
taneously reduce the disk's total mass, outer radius and 
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viscous timescale. These observational facts have gen- 
erally been interpreted as evidence that binaries tighter 
than ~ 100 AU are much less likely to support gas giant 
planet formation. However, follow-up imaging surveys 
have identified some 50 planet-host stars that possess at 
least one stellar companion (e.g., Patience et al. 2002; 
Chauvin et al. 2006; Raghavan et al. 2006; Eggenberger 
et al. 2007; Mugrauer & Neuhauser 2009). In particular, 
it it is worth noting that about 20% of all known plan- 
ets in binary systems have a stellar companion within less 
100 AU, so that planet formation in such an environment 
cannot be considered a rare occurrence. 

In this Letter, I review several key statistical proper- 
ties of PMS and field binary systems that provide insight 
on the planet formation process (Sections 2 and 3). I 
then discuss the implications for the main mechanisms of 
planet formation in binary systems as a function of their 
projected separation (Section 4). In this study, I only 
consider binaries in the 5-1400 AU separation range, for 
which current PMS multiplicity surveys are reasonably 
complete. The tightest binary system known to host a 
planet has a 19 AU separation. Stellar companions be- 
yond 1400 AU are not expected to have much influence 
on planet formation. 

2. PLANET FORMATION IN BINARIES: INITIAL 
CONDITIONS 

2.1. Defining an homogeneous sample of PMS binaries 

In order to draw a broad and homogeneous view of the 
initial conditions for planet formation, I compiled a sam- 
ple of 107 PMS binaries for which deep (sub) millimeter 
continuum observations and/or near- to mid-infrared col- 
ors are available in the literature. The (sub)millimeter 
data are taken from the work of Andrews & Williams 
(2005, 2007); for almost all targets, a la sensitivity of 
15mJy or better at 850/im and/or 1.3mm is achieved. 
The median projected separation in this sample is 92 AU. 
I also defined a comparison sample of 222 PMS stars for 
which no companion has ever been detected. I focus here 
on the Taurus and Ophiuchus star forming regions, the 
only ones for which high-resolution multiplicity, photo- 
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metric and millimeter surveys have a high completeness 
rate. The two clouds contribute an almost equal num- 
ber of binaries to the sample. Furthermore, both regions 
have similar stellar age distributions (median age around 
1 Myr, Ophiuchus being probably slighter younger on av- 
erage than Taurus) and their mass function fully samples 
the 0.1-1.5 Mq range (e.g., Luhman & Rieke 1999; Luh- 
man 2000). Finally, Taurus represents an instance of 
distributed star formation, while Ophiuchus is a more 
clustered environment. These two clouds therefore of- 
fer a global view of the early stages of planet formation 
among solar-type and lower-mass stars. 

2.2. Inner disk presence and survival 

I first address the question of the presence of dust 
in the planet-forming region, namely the innermost few 
AU around each component, within binary systems. To 
probe the presence of an optically thick dusty inner disk, 
I used near- to mid-infrared colors. I selected the follow- 
ing standard thresholds to conclude that a circumstcllar 
disk is present: [3.6] -[8.0] > 0.8 mag, K-N > 1.75 mag, 
K — L > 0.35 mag, a2-i4^m > —1-7 (e.g., Cieza et 
al. 2009; McCabe et al. 2006; Bontemps et al. 2001). 
About 80% of the PMS binaries considered here have 
Spitzer/IRAC colors, which are used whenever available. 

Cieza et al. (2009) have demonstrated that tighter bi- 
naries have a much lower probability of hosting circum- 
stcllar dust. The same effect is observed here in a some- 
what smaller sample. The median separation of bina- 
ries with an inner disk in this sample is about 100 AU, 
whereas that of disk- free binaries is 40 AU. The simplest 
interpretation of this trend is that disks in tight binaries 
are dissipated much faster than in wide systems (Cieza 
et al. 2009, Kraus et al., in prep.). To extend upon this 
previous analysis, I used the two-sided Fischer exact test 
to determine the probability that wide and tight binaries 
have a different proportion of diskless systems, using a 
sliding threshold to split the sample. As shown in Fig- 
ure 1 , the difference is significant at the 2a level or higher 
for a wide range of threshold separations. In particular, 
this analysis reveals that the observed reduced disk life- 
time in binaries only applies to systems that are tighter 
than about 100 A U. On the other hand, there is no statis- 
tical difference between binaries wider than 100 AU and 
single stars. 

2.3. Total mass of the initial dust reservoir 

While near- and mid-infrared emission best traces the 
presence of dust within a few AU of star, only long- 
wavelength flux measurements can probe the total dust 
mass of protoplanetary disks (e.g., Beckwith et al. 1990). 
From the sample defined above, I selected those ob- 
jects which show evidence of an optically thick inner 
disk (as defined above) and have been observed in the 
(sub)millimeter. The median separation in this subsam- 
ple of 44 binaries is 130 AU. While the 850/im survey of 
Ophiuchus is not yet as complete as that of Taurus, the 
existing 1.3mm observations of PMS stars are generally 
less sensitive to cold dust. Since using both wavelengths 
yield similar conclusions but with lower significance for 
the 1.3mm one, I focus here on 850/im measurements. 

As has long been known, tight binaries have a differ- 
ent distribution of submillimetcr fluxes than wide ones, 



with a much lower median flux (13mJy vs 50mJy at 
850/im using a 100 AU separation threshold) and only 
very few high-flux systems (Jensen et al. 1996; Andrews 
& Williams 2005). I compared the distributions of 850/im 
fluxes for tight and wide binaries defined by the same 
sliding threshold as above using the conservative sur- 
vival analysis Peto-Pentrice Generalized Wilcoxon test 
to account for upper limits. I find that wide and tight 
binaries arc different at the 2a level or higher if the sepa- 
ration threshold is in the 75-300 AU range (see Figure 1). 
I therefore conclude that binaries with a projected sepa- 
ration smaller than 300 AU have a substantially reduced 
submillimeter flux. On the other hand, the distribution 
of 850/xm fluxes for wide binaries is indistinguishable 
from that of single stars. 

In past studies, it has been assumed that a reduced 
(sub)millimeter flux necessarily implies a reduced total 
dust mass independently of the disk properties (for in- 
stance, see the prescription used by Andrews & Williams 

2005) . While this is true in general, it is unclear whether 
this assumption is valid for severely truncated disks for 
which optical depth effects may become important. The 
model constructed by Jensen et al. (1996) seems to sup- 
port this hypothesis, but these authors assumed that 
tight binaries are always surrounded by a massive cir- 
cumbinary structure, which we now know is rare. To 
revisit this issue, I have computed a grid of radiative 
transfer models using the MCFOST code (Pinte et al. 

2006) to compute the 850/im flux of a disk with a typ- 
ical S(r) oc r _1 surface density profile, an 0.1 AU inner 
radius and a flaring power law H(r) oc r 1125 . Emission 
from the central star is modeled as a 4000 K, 2L Q pho- 
tosphere and a distance of 140 pc is assumed. The dust 
is assumed to be made of astronomical silicates with a 
fl -3.5 p Gwer j aw s j ze distribution ranging from 0.03/im 
to 1 mm. The only variables in the model are the disk 
outer radius, R ou t, and the total dust mass, Mdust- Fig- 
ure 2 demonstrates that the proportionality between to- 
tal dust mass and submillimeter flux observed for large 
disks breaks down for R out < 30 AU as the disk becomes 
optically thick to its own emission. 

Disk truncation by an outer stellar component is de- 
pendent on the orbital parameters and mass ratio of the 
binary system (Artymowicz & Lubow 1994). It is there- 
fore not possible to uniquely associate a binary separa- 
tion with a tidally-set value of R ou t- The ratio between 
these quantities is typically in the broad 2.5-5 range. 
Systems whose separation is less than 100 AU are there- 
fore expected to possess disks whose outer radius is 40 AU 
or less. In this configuration, total disk masses of at 
least Mj are necessary to produce 850/im fluxes as low as 
~ 20-30 mJy In the sample studied here, about a third 
(6 out of 19) of all binaries that are tighter than 100 AU 
and possess an inner disk have an 850/im flux that is 
higher than 30mJy. Therefore, a significant fraction of 
the circumstellar disks in tight binaries (< 100 AU) are 
massive enough to potentially form gas giant planets, de- 
spite their much lower (sub) millimeter fluxes. 

3. PLANET FORMATION IN BINARIES: END RESULTS 

3.1. Distribution of planetary masses 

Les us turn our attention to mature planetary systems. 
As of this writing, there are 38 exoplanets that are in a 
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system with at least 2 stellar components (using 1400 AU 
as the upper limit for binary separation) , including 5 sys- 
tems with a stellar companions within 25 AU (Mugrauer 
& Neuhauser 2009, and references therein). Most of these 
planet-bearing stars are of solar type. The overall detec- 
tion rates of gas giant planets in binary systems and in 
single stars are undistinguishable (Bonavita & Desidera 
2007). There is marginal evidence that planets in bina- 
ries tighter than ~ 100 AU may be somewhat less fre- 
quent than one would assume based on the frequency of 
planets in wider binaries (by 6.0±2.7%, Eggenberger et 
al. 2008). However, the small sample size, adverse se- 
lection biases and incompleteness of current multiplicity 
surveys are such that it is premature to reach definitive 
conclusions. In any case, we can use this sample to test 
whether the separation of the stellar binary has any in- 
fluence on planet properties. 

Despite an earlier claim for a distinct period-mass dis- 
tribution (Zucker & Mazeh 2002), Bonavita & Desidera 
(2007) have shown that there is essentially no difference 
in the properties of planetary systems around single stars 
and in binary systems. However, a previously unrecog- 
nized trend is evident in Figure 3. While planets covering 
two orders of magnitude in mass can be found in wide 
binaries (as around single stars), systems tighter than 
~ 100 AU appear to host only high-mass, M > Mj, 
planets. To quantify this effect, I used the two-sided 
Fischer Exact test to determine whether close and wide 
binaries (with the usual sliding threshold) have differ- 
ent proportion of high- and low-mass planets. I used 
1.6Mj, the median for all planets known to date 2 , to 
separate low- from high-mass planets. Figure 1 confirms 
that binary systems tighter than about 100 AU produce a 
distribution of planets that is strongly biased towards the 
highest masses. This conclusion is significant at the So- 
level. 

It is important to test whether this trend is not a mere 
consequence of a selection bias, as a close stellar compan- 
ion can alter the detectability of a planet-induced radial 
velocity signal. To evaluate this possibility, I build on 
the "uniform detectability" sample defined by Fischer & 
Valenti (2005) which contains all stars for which close-in 
planets as low mass as 0.3Mj (well below the apparent 
cut-off in mass for planets in tight binaries), as well as 
lMj planets on a 4yr orbit, could be detected. In the 
current sample of binary planet hosts, the proportions 
of stars that belong to the uniform detectability sample 
among binaries tighter and wider than 100 AU are indis- 
tinguishable (5/9 and 24/40, respectively). I therefore 
conclude that the trend discussed above is unlikely to be 
the consequence of a selection bias or of observational 
limitations. 

3.2. Occurrence of planetesimal disks 

An indirect signpost of planet formation is the debris 
disk phenomenon. In these systems, small dust grains are 
produced via the collisions of large solid bodies (Zucker- 
man 2001). Trilling ct al. (2007) observed 69 A- and F- 
type known binaries with Spitzer and found debris disks 
in systems spanning 6 decades in separation. They fur- 
ther suggested that intermediate separation (3-30 AU) 

2 This is also the median for the subsample of all planets within 
multiple stellar systems. 



binaries are substantially less likely to host a debris disks 
than cither tighter or wider systems, although the formal 
significance of this difference is marginal at best. No such 
trend was found by Plavchan et al. (2009), who included 
24 A- through M-type binaries in their own Spitzer sur- 
vey. This latter survey focused on targets that are more 
similar in mass to exoplanet hosts and the PMS popula- 
tion discussed in the previous section. 

I used the two-sided Fischer test to determine whether 
the occurrence of debris disks is indeed different in tight 
and wide binaries, using the same sliding threshold as 
above (see Figure 1). There is no significant difference 
for any value of the threshold in the sample from Trilling 
et al. (2007), nor in a combined sample that also includes 
systems from Plavchan et al. (2009). The combined sam- 
ple contains 52 binaries in the 5-1400 AU range, with a 
median separation of 50 AU, an increase of 15 sources 
from the sole sample of Trilling et al.. In addition, the 
occurrence rates of debris disks in binary systems and 
single stars are very similar (Trilling ct al. 2007). In other 
words, any 0.5-2Mq star, irrespective of the presence of 
a companion (within the 5-1400 A U range studied here), 
may experience the early phases of planet formation up 
to the planetesimal stage. 

4. DISCUSSION AND IMPLICATIONS 

This analysis has revealed a clear dichotomy between 
tight and wide binaries. Systems with separation > 
100 AU are indistinguishable from single stars as far as 
the initial conditions and end product of planet forma- 
tion are concerned. The only caveat to this statement 
is the possibility of mild disk truncation in 100-300 AU 
systems, but most disks in these systems retain a mass 
reservoir that is sufficient to build up gas giant planets. 
On the other hand, planet formation in binaries with 
separations < 100 AU is characterized by a much shorter 
clearing timescale for the protoplanetary disks and a 
strong bias towards high-mass planets. Despite these 
differences, planetesimals and mature planetary systems 
appear to form at roughly the same frequency as around 
other stars. Furthermore, while protoplanetary disks are 
more compact in tight systems because of truncation, 
a significant fraction of them possess large mass reser- 
voirs (at least several times Mj). Taken together, these 
results suggest that planet formation in binaries tighter 
than 100 A U proceeds through a different, but not much 
less frequent, mechanism compared to wide binaries and 
single stars. 

The shorter disk lifetime in tight binaries makes it ex- 
tremely difficult to form gas giant planets through the 
core accretion model, especially if the final planets are 
particularly massive. Rather, this combination of ob- 
served trends supports an abrupt process to form planets 
in tight binaries, such as the disk fragmentation model. 
Indeed, this mechanism can be extremely efficient in the 
case of a compact, massive protoplanetary disk which 
is naturally prone to gravitational instability. Further- 
more, gravitational perturbations induced by a close stel- 
lar companion can trigger the instability even though 
the disk itself is not unstable to its own gravity (Boss 
2006). On the other hand, considering the long survival 
timescale and slim chances of gravitational instabilities, 
disks located within wide binaries and around single stars 
are good candidates to form planets via the core accre- 
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tion model in their inner regions (e.g., Boley 2009). 

While a violent process is most likely responsible for 
the formation of planets in tight binaries, it is however 
unclear whether all planets in wide binaries form through 
a single mechanism. Indeed, it is also conceivable that 
high-mass planets (> Mj) mostly form via disk frag- 
mentation, while lower mass planets are preferentially 
the result of core accretion. This scenario would natu- 
rally alleviate the difficulty of the core accretion model 
to form the highest-mass planets in less than a few Myr. 
This hypothesis has the additional advantage that it 
could also apply to tight binaries. Indeed, since a stel- 
lar companion located within less than 100 AU dramat- 
ically shortens the disk lifetime, core accretion is essen- 
tially prevented from occurring, accounting for the ab- 
sence of low-mass (< 2Mj) gas giant planets in tight 
binaries. Planetesimals can presumably form in either 
scenario, accounting for the observations regarding the 
debris disks phenomenon. In summary, it remains to be 
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determined whether the trends discussed here indicate 
an actual dichotomy between the main planet formation 
theories or a mere change of the relative importance of 
the two models as a function of the location of the stellar 
companion. Improving the statistical significance of the 
various trends discussed here and determining the exact 
properties of disks within tight PMS binaries will help 
shed further light these two possibilities. 
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TABLE 1 

Sample of Taurus and Ophiuchus protoplanetary disks considered in this study. 



Target 


Alt. Name Sep. 


IR color a 


Inner disk? 




Pl.3mm 


References 15 




[AU] 


[mag] 




[mjy] 


[mjy] 




Taurus-Auriga 


CZ Tau 


46 


3.41 


Y 


< 9 


< 30 


1,2 


DD Tau 


79 


1.99 


Y 


< 42 


17 


1,2 


DF Tau 


13 


1.37 


Y 


8.8 


< 25 


3,2 


DI Tau 


17 


0.76 


N 






3 


DK Tau 


350 


1.27 


Y 


80 


35 


3,2 



a Unless otherwise indicated by a note, the infrared color is the Spitzer/IRAC [3.6] — [8.0] color. 

b References: 1) Luhman et al. (2006); 2) Andrews & Williams (2005); 3) Hartmann et al. (2005); 4) Cieza et al. (2009); 5) McCabe et al. 
(2006); 6) Guilloteau et al. (1999); 7) Kenyon & Hartmann (1995); 8) Duchene et al. (2010); 9) Jensen & Akeson (2003); 10) Jensen et al. 
(1996); 11) Andrews & Williams (2007); 12) Bontemps et al. (2001). 
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100 

Tight/wide binary separation threshold (AU) 

Fig. 1. — Probability of several null hypotheses that compare the observed properties of tight and wide binaries as a function of the radius 
used to separate both categories. The following null hypotheses are tested: planets in tight and wide binaries have the same frequency 
of high mass planets, namely above the median for all known extrasolar planets (red solid curve); tight and wide binaries have the same 
probability of hosting optically thick dust within a few AUs of the central stars (blue dashed curve); the distributions of 850/im fluxes for 
tight and wide binaries with evidence of optically thick inner disks are indistinguishable (green dot-dashed curve); the frequency of the 
debris disk phenomenon is the same for tight and wide binaries (purple triple- dot- dashed curve). The horizontal dotted lines represent the 2 
and 3cr confidence levels (top and bottom, respectively). Notwithstanding small number statistics (responsible for the jagged appearance of 
the curves), the first three null hypotheses can be rejected at the 99.5% confidence level. Only the debris disk phenomenon is independent 
of the binary separation. 
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Fig. 2. — Model 850^im flux (in mjy) for a typical flared disk (see model parameters in section 2.3) with varying total dust mass M^ust 
and outer radius Rout, viewed at an inclination of 60° (see Section 2.3 for more detail). The horizontal dotted line represent a total disk 
mass Md = Mj, assuming a 100:1 gas-to-dust ratio. Compact disks, with R ou t < 10 AU must be massive to yield -FssO/im > 20mjy. 
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Fig. 3. — Distribution of masses (Msini) for planets in binary systems as a function of the binary separation. Note the deficit of planets 
with Afsini < Mj for binaries tighter than 100 AU. The horizontal dotted line represent the median mass for all known extrasolar planets 
to date (including among single stars). 



